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The design of nano-antennas is so far mainly inspired by radio-frequency technology. However, 
material properties and experimental settings need to be reconsidered at optical frequencies, which 
entails the need for alternative optimal antenna designs. Here a checkerboard-type, initially random 
array of gold cubes is subjected to evolutionary optimization. To illustrate the power of the approach 
we demonstrate that by optimizing the near-field intensity enhancement the evolutionary algorithm 
finds a new antenna geometry, essentially a split-ring/two- wire antenna hybrid which surpasses by 
far the performance of a conventional gap antenna by shifting the n=l split-ring resonance into the 
optical regime. 

PACS numbers: 84.40.Ba, 73.20.Mf, 02.60.-x, 78.67.Bf 



Light-matter interaction, i.e. absorption and emission 
of light as well as the control of its spectral and direc- 
tional properties, can be optimized by means of antenna- 
like plasmonic nano structures [1,2]. This is of immediate 
importance in diverse fields of research ranging from so- 
lar energy conversion [3], photocatalytic [4] and sensing 
applications [5] to single-particle manipulation [6, 7] and 
spectroscopy [8] as well as quantum optics and commu- 
nication [9-12]. 

RF- antenna designs are usually optimized for thin, in- 
finitely good conducting wires that only support surface 
currents and are typically fed by transmission lines con- 
nected by infinitely narrow gaps [13]. For antennas at op- 
tical frequencies the general operation conditions deviate 
substantially from such ideal behaviour: (i) Antenna wire 
diameters are comparable to the electromagnetic pene- 
tration depth into the wire material leading to volume 
currents [14]. In the case of noble metals, such wires 
therefore exhibit plasmon resonances in the visible spec- 
tral range resulting in a reduced eflFective wavelength of 
wire waves [15]. (ii) Feeding (excitation) of optical anten- 
nas is often achieved by focused laser beams or quantum 
emitters, (iii) high-frequency-related eflFects such as the 
'kinetic inductance' become significant [16]. It can there- 
fore not be taken for granted that RF-inspired antenna 
designs, like dipole [17], bow tie [18, 19] and Yagi-Uda 
antennas [20, 21], represent 'optimal' geometries also at 
optical frequencies, although they provide a reasonable 
performance. 

Evolutionary algorithms (EAs) find optimized solu- 
tions to highly complex non-analytic problems by creat- 
ing subsequent generations of individuals coded by their 
respective genomes that compete for the right to pass 
on their properties, according to a fitness parameter [22]. 
These optimized solutions can then be analyzed to foster 
the understanding of underlying physical principles. Evo- 
lutionary optimization has successfully been applied in 
various fields of research, including pulse shape optimiza- 
tion in coherent control of chemical reactions [23] and 



field localization in plasmonic structures [24, 25]. Fur- 
thermore, evolutionary optimization has been used to aid 
the development of radio- wave antennas [26, 27]. First 
attempts to employ such methods to find improved plas- 
monic nanostructures have also been undertaken [28-31], 
however, the investigated configurational space remained 
very limited. 

Here we employ the method of evolutionary optimiza- 
tion in a more general setup to find improved plas- 
monic antenna structures that, in terms of near-field 
intensity enhancement (fitness parameter), outperform 
the best radio- wave- type reference antennas by a fac- 
tor of two. Analysis of the fittest antenna reveals that 
it is a split-ring/two- wire antenna hybrid (Split-Ring- 
Antenna), which merges features of the fundamental 
magnetic resonance of a split ring with the fundamen- 
tal electric resonance of a linear dipole antenna both in 
the visible wavelength regime. 

Methods As fitness parameter we choose the normal- 
ized near-field intensity enhancement in the focus of an 
illuminating Gaussian beam (A^ =647 nm, NA= 1, 4.3 fs 
pulse duration, 144 nm bandwidth). A genetic represen- 
tation of complex-shaped thin-film nano antennas is real- 
ized by composing structures (matrix antennas) from dis- 
crete gold cubes with fixed dimensions (10x10x11 nm^) 
positioned on a 21x21 square matrix in vacuum oriented 
perpendicular to and centered on the optical axis of the 
Gaussian beam. The fields of such structures can be de- 
scribed by local Maxwell equations. 

The resulting configuration space of about 4-10^^^ dif- 
ferent individual structures ensures geometrical variety 
but is impossible to explore by brute force methods, since 
the evaluation of an individual structure takes about 20 
minutes. The size of the focal spot and the area occupied 
by the gold cube matrix are comparable (see Fig. SI [32]). 
An example showing a bow-tie antenna represented in a 
5x5 array is depicted in Fig. 1. The genetic information 
is represented in a unique binary code, where Matrix ele- 
ments are set to '1', if occupied by a gold cube and to '0' 



2 



when empty (for a calculation of the number of physical 
redundant structures see Fig. S2, S3 and related discus- 
sion) . 

We solve Maxwell's equations using the finite- 
difference time-domain (FDTD) method [33] (FDTD So- 
lution, Lumerical Inc., Canada). Due to the discretiza- 
tion in space (Yee Cells, dx = 1 nm), adjacent gold cubes 
are fully connected to each other. We also obtain con- 
ductive bridges between neighbouring gold blocks that 
contact each other via their edges. All simulation objects 
were shifted by half a Yee cell size in x- and ^/-direction to 
ensure identical dimensions of gold cubes and voids. In 
order to describe the dielectric function of gold in a suf- 
ficiently wide spectral range we use the data by Johnson 
and Christy fitted by an analytical model [34]. 
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FIG. 1. (color online) Genetic representation of a matrix an- 
tenna. Left: artistic 3D-view of an example bow-tie nano an- 
tenna consisting of 8 gold cubes. Right: top view indicating 
the transition to a 5x5-matrix representation, where each '1' 
denotes the presence of a 10x10x11 nm^ block of gold. 



The evolutionary algorithm is implemented in MatLab 
[35]. It uses generations consisting of 20 or 30 individ- 
ual matrix antennas. The five best structures according 
to the fitness parameter are selected as parents for the 
next generation. With carefully chosen mechanisms for 
crossover and mutation of the genomes, consecutive gen- 
erations constantly improve in the sense that their fitness 
parameter increases. To create descendants in a first step 
one of the five parents is selected by roulette wheel selec- 
tion [22] with a probability that is proportional to its fit- 
ness (parent A). Three methods, i.e. creation of random 
structures, mutation as well as linear and spiral genome 
crossing, are then applied until a total of 20 (or 30) new 
individuals have been generated for the next generation 
(see Fig. S4). 

After a sufficient amount of generations has been simu- 
lated, a so called toggle plot analysis is performed, which 
consists of running 21x21 simulations in which every 
block is toggled individually. Colour-coding the block po- 
sitions according to the magnitude of the fitness changes 
associated with the individual toggle event shows the rel- 
ative importance of single blocks, the potential for further 
improvement of a matrix antenna, and also eventually 
produces new individuals with enhanced fitness [36]. 



Results The wavelength coincides with the reso- 
nance of a linear dipole nano antenna consisting of two 
end-to-end aligned 46x30x11 nm^ gold rods (width = 
3 cubes, height = 1 cube) separated by a 10 nm gap 
(Fig. 2(a), top panel) which serves as a reference struc- 
ture. It exhibits a resonant normalized near-field inten- 
sity enhancement of about 1800 in the center of its feed 
gap. Other geometries, such as bow-tie antennas, were 
also tested but do not yield higher fitness. 

In the following we discuss the fittest structure ob- 
tained by running the EA for 100 generations with 20 in- 
dividuals each (see Fig. S5 for development of geometry 
and fitness parameter), a subsequent toggle plot analy- 
sis and further 30 generations with 30 individuals each, 
starting with combinations of the best five structures ob- 
tained from the toggle plot analysis. 

The best matrix antenna structure (Fig. 2 (a), lower 
panel) exhibits a remarkably high fitness, as indicated 
by its near- field spectrum in Fig. 2(b) which is recorded 
in the optimization point (indicated as blue dot in (a)) 
after a broadband excitation. Its maximal near-field in- 
tensity enhancement of 3500 is nearly twice as high as 
that of the reference antenna. Both spectra show single 
nearly Lorentzian peaks (Q=20 and 23, respectively; see 
supplementary information for a further discussion of the 
Q-factor). 

According to the reciprocity theorem [1, 2] the opti- 
mized antenna should also improve the radiative proper- 
ties of a quantum emitter positioned in the spot of highest 
field enhancement. Indeed for the reference antenna we 
find a radiation enhancement of 2126 and a radiation effi- 
ciency of 0.255, while for the fittest antenna the radiation 
enhancement is 4271 with a radiation efficiency of 0.268. 
Surprisingly, the directivity of the fittest antenna remains 
very similar to that of the reference antenna despite its 
complex shape (see Fig. S6 and related discussion). 

The fittest matrix antenna exhibits three noticeable 
geometrical features: 

(i) a small gap in the center between two compact rod(- 
like) structures, being slightly displaced in y-direction 
with respect to the observation point. 

(ii) a single gold block directly below the gap which 
creates a current path connecting the rod-like structures 
and 

(iii) a seemingly random arrangement of gold blocks 
further away from the center. 

It is important to note that the optimal structure found 
by a genetic algorithm depends on both the available 
primitive elements as well as on possible boundary con- 
ditions. Using a different block size or imposing bound- 
ary conditions, e.g. limiting the gap width, would lead 
to different fittest structures. 

We now consider the near-field intensity enhancement 
maps of both reference and matrix antenna in Fig. 2(c). 
The small displacement of the rod-like structures in- 
creases the near-field intensity enhancement by a small 
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FIG. 2. (color online) Comparison at at A =647 nm between a resonant linear dipole nano antenna build from two 10 nm 
separated rectangular arms of 46x30x11 nm^ and a structure obtained with the evolutionary algorithm described in the text, (a) 
shows the geometry of both structures from the -\-z direction. The blue spot denotes the position of the near field optimization 
by the EA. The spectra in (b) are taken at this marked position during a broadband gaussian excitation, (c) shows the 
logarithmic near field intensities at A =647 nm when the structures are illuminated by a monochromatic gaussian focus with 
NA = 1. The scales are normalized and not comparable. In (d) the the change of near field enhancement at the optimization 
position is shown for each single block, when it is toggled, (e) is a zoom of the centre part of the EA antenna, showing strength 
and direction of the currents. 



factor because of the proximity of the corners of the rod- 
like structures to the point of optimization. However, 
this alone by far cannot explain the observed increase 
of the near- field intensity enhancement. The achievable 
enhancement by displacing the reference antenna in a 
similar way amounts to a factor of 1.1. 

The result of a toggle plot analysis is displayed in 
Fig. 2(d). It indicates that changing individual blocks 
does not yield considerable additional near-field inten- 
sity enhancement, but rather a severe reduction. We 
therefore believe that the structure's fitness is close to 
a (local) maximum in the configuration space. The by 
far strongest reduction of fitness occurs when toggling 
gold blocks near the center. This indicates that the com- 
pact structure in the proximity of the gap is dominating 
the field enhancement and is most critical for achieving 
the observed performance. Assuming that it is excited at 
an Eigenmode, this also explains the appearance of a sin- 
gle narrow Lorentzian resonance. As apparent from the 
toggle plot, the random structures far away from the cen- 
ter do hardly influence the field enhancement in the gap. 
Nevertheless, it is possible that collective effects of the 
peripheral blocks do influence the fitness of the structure 
to some extend which due to the inherent complexity will 
not be further discussed in detail. 

Of particular interest is the single gold block below the 
gap. It provides a current path via the cube edges be- 
tween the two rod-like structures that form the gap as 
can be seen by taking a closer look at the currents in the 



central part of the matrix antenna in Fig. 2(e). Surpris- 
ingly, we find that removing this block severely lowers the 
fitness of the resulting matrix antenna instead of increas- 
ing it. Closer inspection of Fig. 2(e) reveals two partic- 
ular current paths, one located in the rod-like structures 
corresponding to a bonding linear dipolar nano antenna 
mode, but also a second one, which flows from one up- 
per gap edge through the connecting gold block to the 
other upper gap edge, corresponding to a fundamental 
split ring mode. 

In order to better understand the effects that lead to 
the increased near-field intensity enhancement, in the fol- 
lowing we study a reduced model system, i.e. a mixture 
of a split-ring and a linear two- wire antenna called split- 
ring antenna (SRA) that retains the important features 
of the fittest matrix antenna but can be described by a 
small number of freely tunable parameters. Its geome- 
try is depicted in Fig. 3(a). The structure can be inter- 
preted either as linear two-wire antenna with an asym- 
metric short circuit, a split-ring resonator with attached 
wires or a long single nano wire that is deformed in a 
particular way. 

The near-field intensity enhancement of the SRA in the 
center plane is depicted in Fig. 3(b), showing a strong 
field concentration towards the open side of the SRA gap 
as it was already observed in the best matrix antenna. 
Fig. 3(c) compares the spectra of the resonances of a 
SRA with a non short-circuited dipolar antenna of iden- 
tical arm cross section, gap size and resonance frequency. 
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FIG. 3. (color online) (a) Geometry and dimensions of the examined split-ring antenna, (b) Near field intensity of the same 
structure, showing the very concentrated fields in the open end of the gap. The blue point denotes the the point of measurement 
for the spectrum figured in (c). The Q ^ 25 for the SRA (red) is comparable to the reference antenna (black), (d) is an overlay 
of the current intensities and direction inside the material, and the magnetic near fields outside the structure. 



Both spectra were obtained at the point of highest near 
field enhancement along the y-ax.is. Also in the present 
model system the split-ring antenna surpasses the classi- 
cal dipole antenna design in terms of maximum near-field 
intensity enhancement by a factor of 2. 

Also in the model system the current pattern can be 
decomposed into a fundamental (n=l) split-ring mode 
[37] and a dipolar current in each antenna arm which 
is running 180° out-of-phase to the current in the short 
circuit (Fig. 3 (d)), adding to the charge accumulation in 
the upper part of the gap and thus increasing the near- 
field intensity enhancement. Since the resonance is in 
the optical regime, the SRA is a way to circumvent the 
limitation of pure split-ring resonances to wavelengths 
above 900 nm due to the kinetic inductance [16]. The 
SRA represents a magnetic dipole in the visible, showing 
magnetic fields (Fig. 3 (d)), which are only by a factor of 
2.5 weaker than those of the isolated split-ring resonator 
at its resonance wavelength of 908 nm. 

To confirm the shift of the fundamental split-ring mode 
from the infrared into the visible spectral range, we place 
two gold bars (35x30x11 nm^) separated by a 5 nm non 
conductive gap from the plain split-rings ends. The ex- 
tinction cross section is examined while connecting the 
gold bars with the split-ring via successively thicker gold 
bridges at the center of the gap (see Fig. 4 and inset). 
In the unconnected geometry the fundamental split-ring 
resonance is shifted into the red to about 1200 nm due to 
the capacitive coupling across the gap. As the connection 
grows thicker (increasing (i), the fundamental split-ring 
mode is shifted by more than 500 nm from the infrared 
into the visible. The shift of the resonance is of similar 
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FIG. 4. Wavelength shift of the extinction cross section res- 
onance for the combined two-wire antenna and the n — 1 
split ring resonance, as a conductive link between two an- 
tenna rods and a split ring (see inset) increases in thickness 
(marked thick dashed). Also visible are the bonding and an- 
tibonding n = 3-modes, which also shift slightly into the blue 
and get very weak (marked thin dashed). The sketches at 
the bottom show the geometries of the fully disconnected and 
fully connected structures together with the position and sign 
of their mode's charge density maxima. 



origin as the emergence of a charge transfer mode for a 
dipolar antenna with a conductive bridge [38]. However, 
here the fundamental split-ring mode is not disappearing 
but its phase is inverted (compare the charge distribu- 
tions sketched in the lower panel of Fig. 4). For very 
thin conductive bridges both modes exist simultaneously 
(see also [38]) and cancel out each other, leading to a dip 
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in the extinction cross section. 

Conclusion We show that by using the method of evo- 
lutionary optimization in a large parameter space, high- 
fitness plasmonic antennas can be found within a reason- 
able amount of time. The method can be adapted to a 
large variety of fitness parameters in order to optimize 
plasmonic structures for various purposes. Besides di- 
rectly yielding optimized structures, careful analysis of 
the working principles of the resulting geometries may 
provide new design strategies for high-performance plas- 
monic nano structures. 

In the present example we obtain an increase of the 
near-field intensity by nearly a factor of two caused by the 
intriguing cooperation of a fundamental split-ring mode 
and dipole antenna resonance. We have used the latter 
principle to devise a novel antenna design which addition- 
ally exhibits very large magnetic fields at optical frequen- 
cies. In particular, the fundamental split-ring resonance 
is shifted into the visible spectral range because of the 
formation of a charge-transfer-like hybrid resonance with 
two rods of a dipolar antenna. The method can be fur- 
ther adapted to include geometrical constraints imposed 
by micro fabrication and therefore can lead to structures 
that can directly be implemented in practical applica- 
tions. 

Acknowledgements Markus Kiunke wants to thank 
Christoph Br lining, Thorsten Feichtner and Oleg Selig 
want to thank Jord C. Prangsma and Paolo Biagioni for 
fruitful discussions. Financial support by the DFG is 
gratefully acknowledged (HE5618/1-1). 



[1] P. Bharadwaj, B. Deutsch, and L. Novotny, Adv. Opt. 

Phot. 1, 438-483 (2009). 
[2] P. Biagioni, J.-S. Huang, and B. Hecht, Rep. Prog. Phys. 

75, 024402 (2012). 
[3] H. A. Atwater and A. Polman, Nat. Mater. 9, 205 (2010). 
[4] Z. Liu, W. Hou, P. Pavaskar, M. Aykol, and 

S. B. Cronin, Nano Letters 11, 1111 (2011), 

http://pubs.acs.org/doi/pdf/10.1021/nll04005n. 
[5] J. Becker, A. Triigler, A. Jakab, U. Hohenester, and 

C. Sonnichsen, Plasmonics 5, 161 (2010). 
[6] M. L. Juan, R. Gordon, Y. Pang, F. Eftekhari, and 

R. Quidant, Nat. Phys. 5, 915 (2009). 
[7] W. Zhang, L. Huang, C. Santschi, and O. J. F. Martin, 

Nano Letters 10, 1006 (2010). 
[8] A. Kinkhabwala, Z. Yu, S. Fan, Y. Avlasevich, 

K. Muellen, and W. E. Moerner, Nat. Photon. 3, 654 

(2009). 

[9] a. V. Akimov, a. Mukherjee, C. L. Yu, D. E. Chang, 

a. S. Zibrov, P. R. Hemmer, H. Park, and M. D. Lukin, 

Nature 450, 402 (2007). 
[10] R. Kolesov, B. Grotz, G. Balasubramanian, R. J. Stohr, 

A. A. L. Nicolet, P. R. Hemmer, F. Jelezko, and 

J. Wrachtrup, Nat. Phys. 5, 470 (2009). 
[11] J.-S. Huang, T. Feichtner, P. Biagioni, and B. Hecht, 

Nano Letters 9, 1897 (2009). 



[12] Z. Jacob and V. M. Shalaev, Science 334, 463 (2011). 

[13] C. A. Balanis, Proc. IEEE 80, 7 (1992). 

[14] J. Dorfmueller, R. Vogelgesang, W. Khunsin, C. Rock- 

stuhl, C. Etrich, and K. Kern, Nano Letters 10, 3596 

(2010). 

[15] L. Novotny, Phys. Rev. Lett. 98, 266802 (2007). 

[16] J. Zhou, T. Koschny, M. Kafesaki, E. N. Economou, 

J. B. Pendry, and C. M. Soukoulis, Phys. Rev. Lett. 

95, 223902 (2005). 
[17] P. Miihlschlegel, H.-J. Eisler, O. J. F. Martin, B. Hecht, 

and D. W. Pohl, Science 308, 1607 (2005). 
[18] P. Schuck, D. Fromm, A. Sundaramurthy, G. Kino, and 

W. Moerner, Phys. Rev. Lett. 94, 017402 (2005). 
[19] J. N. Farahani, D. W. Pohl, H.-J. Eisler, and B. Hecht, 

Phys. Rev. Lett. 95, 017402 (2005). 
[20] T. H. Taminiau, F. D. Stefani, and N. F. Van Hulst, 

Opt. Expr. 16, 10858 (2008). 
[21] A. G. Curto, G. Volpe, T. H. Taminiau, M. P. Kreuzer, 

R. Quidant, and N. F. van Hulst, Science 329, 930 

(2010). 

[22] S. Sivanandram and S. Deepa, Introduction into genetic 
algorithms (Springer Verlag Berlin Heidelberg, 2008). 

[23] T. Baumert, T. Brixner, V. Seyfried, M. Strehle, and 
G. Gerber, Appl. Phys. B: Lasers and Optics 65, 779 
(1997). 

[24] M. Aeschlimann, M. Bauer, D. Bayer, T. Brixner, F. J. 
Garcia de Abajo, W. PfeifFer, M. Rohmer, C. Spindler, 
and F. Steeb, Nature 446, 301 (2007). 

[25] M. Aeschlimann, M. Bauer, D. Bayer, T. Brixner, 
S. Cunovic, F. Dimler, A. Fischer, W. Pfeiffer, 
M. Rohmer, C. Schneider, F. Steeb, C. Striiber, and 
D. V. Voronine, Proc. Natl. Acad, of Sci. USA 107, 5329 
(2010). 

[26] H. Huang, A. Hoorfar, and S. Lakhani, Antennas 

and Propagation Society International Symposium, 2007 

IEEE 1, 1609 (2007). 
[27] M. F. Pantoja, A. R. Bretones, S. Member, R. G. Martin, 

and S. Member, IEEE Trans. Antennas Propag. 55, 1111 

(2007). 

[28] P. Ginzburg, N. Berkovitch, A. Nevet, I. Shor, and 
M. Orenstein, Nano Letters 11, 2329 (2011). 

[29] S. Kessentini, D. Barchiesi, T. Grosges, and M. de la 
Chapelle, Evolutionary Computation (CEC), 2011 IEEE 
Congress on 1, 2315 (2011). 

[30] C. Forestiere, M. Donelli, G. F. Walsh, E. Zeni, G. Miano, 
and L. D. Negro, Opt. Lett. 35, 133 (2010). 

[31] C. Forestiere, A. J. Pasquale, A. Capretti, G. Mi- 
ano, A. Tamburrino, S. Y. Lee, B. M. Reinhard, and 
L. Dal Negro, Nano Letters 12, 2037 (2012). 

[32] See Supplemental Material at [URL will be inserted by 
publisher] . 

[33] A. Taflove and S. C. Hagness, Computational electro- 
dynamics: the finite- difference time-domain method, 3rd 
ed., edited by A. Taflove (Artech House, Inc., 2005). 

[34] P. G. Etchegoin, E. C. L. Ru, and M. Meyer, J. Chem. 
Phys 125, 164705 (2006). 

[35] All necessary code written in Lumerical script language 
or MatLab can be made available upon request. 

[36] This is the brute force realization of the hill climber al- 
gorithm on our setup, neglecting pair- and higher-order 
correlations. 

[37] C. Rockstuhl, F. Lederer, C. Etrich, T. Zentgraf, J. Kuhl, 
and H. Giessen, Opt. Expr. 14, 8827 (2006). 



6 



[38] M. Schnell, A. Garcia-Etxarri, A. J. Huber, K. Crozier, 
J. Aizpurua, and R. Hillenbrand, Nat. Photon. 3, 287 
(2009). 



Evolutionary optimization of plasmonic nanostructures 
— Supplementary material — 

Feichtner Thorsten, Oleg Selig, Markus Kiunke, and Bert Hecht 
Nano- Optics & Biophotonics Group, Department of Experimental Physics 5, 
Rontgen Research Center for Complex Material Research (RCCM), 
Physics Institute, University of Wiirzburg, 
Am Hubland, D-97074 Wiirzburg, Germany 



SIZE RELATION OF MATRIX ANTENNA AND FOCAL AREA 

The size of the matrix and the gold cubes were chosen such that the whole structure is 
illuminated by the excitation focus (Fig. SI) and every cube is subjected to about the same 
light intensity. 




FIG. 1. Extension of the 21x21-matrix compared to the field profile of the 647 nm gaussian focus 
with NA 1. 
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NUMBER OF NON REDUNDANT BINARY ODD-SIZED n x n-MATRICES 
WITH ONE SYMMETRIC AND ONE ANTISYMMETRIC AXIS IMPOSED BY 

THE EXCITATION POLARIZATION 

We consider a binary n x n-matrix A with odd n, A^^y G [0, 1] and A(^+i)/2,(n+i)/2 = 0. It 
describes a plasmonic nano structure which is oriented in the x-y-plane centered on the z- 
axis and thus to the optical axis of the system. The hnear excitation polarization is oriented 
along the x-axis resulting in electric fields, that are symmetric with respect to the x-z- 
plane and antisymmetric with respect to the ^-z-plane. Consequently, structures that diflFer 
among each others only by a respective symmetry operation will result in the same physical 
behaviour, although the matrix apparently changes. For matrix sizes with n ^ oo one can 
estimate that about 3/4 of the possible structures are physically redundant. Considerable 
computational eflFort can therefore be avoided by checking each newly generated structure 
for redundancy with already simulated structures. 

First we calculate the number of structures with special symmetries within the set of all 
matrices as defined above. 
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a) axial symmetry b) point symmetry c) axial & rotation 

symmetry 

FIG. 2. At the example of 5 x 5 matrices, the necessary areas to describe structures with given 
symmetry axes (grey) are marked in red. 



Axial symmetry 

In the case of axial symmetry the values in one row containing the center (n — 1 elements) 
and a half matrix (rectangle of n • (n — l)/2 elements; see fig. 2a) ) completely define a 
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structure. The overall number Na of independent structures is thus: 

Na = 2('^-i)+"-'^ = 2("-'Ht+i) 

Point symmetry 

In the case of point symmetry, the values in half a row ((n — l)/2 elements) and half a 
matrix completely describe the structures (see fig. 2b) ), yielding 

n-l I n(n-l) , n^-l 

A^P = 2^+^^ = 2^ (2) 

different configurations. 

Point and axial symmetry 

Every matrix, which is symmetric to both axis is also point symmetric, and every matrix 
which is symmetric to the origin and to one axis, is also symmetric to the other axis. Thus 
only one overlapping set of all three single symmetry sets has to be computed. It is defined 
by a quarter matrix including the on-axis elements and excluding the center block (see Fig. 
SI (c)): 

Nap = 2m'-' (3) 
Counting 

Each matrix, which does not have one of the described symmetries can be transformed 
into three physical identical matrices by mirroring at the x- and ^/-axis. Thus the number of 
non physically redundant non symmetric structures is their complete number N^s divided 
by four. N^s can be computed as the number of all possible structures A^oo = 2^^~^ minus all 
symmetric structures. As there are two mirror axes, the number of axis symmetric matrices 
Na has to be subtracted twice and the number of point symmetric structures once. The 
set of both point- and axis-symmetric matrices has then been subtracted three times and is 
re-added 2 times therefore. This leads to: 

A^oo - 2A^^ -Np + 2Nap 
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Each axis and point symmetric structure can be mirrored into exactly one physical redun- 
dant but geometrical different structure. So to get the number of physical unique matrices, 
their number has to be divided by two after subtracting the structures with both symmetries. 



N, 



2- Na + Np-3- Nap 



(5) 



The structures with axes and point symmetry have to be counted completely, since no 
symmetry operation changes them neither geometrically nor physically. 

The resulting complete number N of physically unique square matrices with odd side 
length n is: 



N = Nns + Na,p + Nap 



(6) 
(7) 



An evaluation of this formula for increasing matrix sizes is depicted in Fig. S3. It 
shows, that for matrix sizes of 5 and bigger 75 percent of all structures are redundant. A 
considerable speed enhancement of the evolutionary algorithm is achieved by checking newly 
generated matrices for redundancies with already evaluated structures. 
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FIG. 3. Amount of non redundant structures for the setup described in the text. For larger matrix 
sizes the limit of the fraction of not physically redundant square matrices is 0.25 (blue dashed line). 



SPEED AND CONVERGENCE 

To reduce the computational effort, simulations were terminated after 35 fs, when more 
than 98% of the initial power has left the simulation volume. Although such a short Simula- 
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tion is not sufficient to reach a very high absolute accuracy, it is stiU possible to determine 
the relative hierarchy of simulated antennas in terms of fitness parameters such as near-field 
intensity. Considering that a full FDTD simulation of an antenna structure generally takes 
70 - 80 fs to fully converge, the overall simulation time is cut by about a factor of two. The 
best individuals are re-simulated until full convergence in order to obtain accurate results. 



EVOLUTIONARY ALGORITHM 



Methods of creating new individuals 



The EA uses different mechanisms to generate individuals for generation n + 1 after 
evaluation of generation n: 



a) 




b) 



c) 



paO.l 








FIG. 4. Illustration of methods used for creating descendants from eligible parents, (a) Mutation 
scheme. Each block is toggled with a certain probability, (b) and (c): Crossing of genomes after 
(b) column-by-column separation as well as (c) spiral-type separation. 



(i) Random: Completely random structures are generated with a 50% chance of each 
block being gold or void. These structures are independent of the parent pool and intended 
to introduce diversity into the genetic material. 
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(ii) Mutation: Each block of parent A is toggled between gold and void with a chance 
of 10% (S4 (a)). This allows good structures to enhance further without intermixing with 
other genomes. 

(ill) Linear and spiral crossing: The line-by-line (spiral-) encoded genomes [see S4 (b) 
and (c)] of two parents, A and are crossed by combining the first half of the genome of 
A with the second half of the genome of B preserving the overall genome length [1]. The 
point at which the individual genomes are split is chosen randomly. This method allows the 
combination of the left (inner) part of parent antenna A with the right (outer) part of a 
second parent B [reddish areas in S4 (b) and (c)]. 

Each method to generate a subsequent individual is applied with equal probability. 
Method (iii) requires a second parent B that also is chosen from the pool of the remaining 
4 fittest individuals in the same way as A. The new individuals are checked if they or a 
physically redundant structure already has been simulated and replaced in that case with 
newly generated individuals (compare with SI). 

Development of higher-fitness individuals by means of evolutionary optimization 

Figure S5 displays the progress of the evolutionary optimization as a function of the 
number of generations. Fig. S5 (a) shows the near- field intensity enhancement in the gap of 
all 2000 simulated individuals of the first run. It is important to note, that the individuals 
have only been simulated for 35 fs and therefore their fitness parameter is underestimated 
typically by 10-20%. The first generation contained 20 random individuals with field en- 
hancements between 0.2 and 43.8. Figure 5 (b) shows, to the very left, the probability 
of each matrix space of the best 10 structures of the first generation to be occupied with 
gold (Gen: 1). The low saturation reflects the randomness of each single matrix position. 
While the algorithm develops, the near-fleld intensity enhancement gradually increases and 
surpasses the enhancement of the reference antenna in generation 42. 

The probability plots in Fig. 5 (b) clearly illustrate the development of some structural 
features: left of the gap a small arm with 40x20 nm size develops, while on the right side 
a more triangular structure emerges. On the way to a better structure the individuals in a 
single generation get more similar within the flrst 10 generations and show only small dif- 
ferences between the best single structures henceforth. The development of some important 
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a) 
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# of simulation 

Gen: 1 Gen: 34 Gen: 55 Gen: 88 




U (always empty) Probability per generation (aiwaysgoid) 1 



FIG. 5. Overview over the evolution process, a) shows the development of the near- field intensity 
enhancement as the number of generations progresses. Twenty simulations make up one generation. 
Number 1, 34, 55, and 88 are denoted by grey bars for further discussion in b). The red curve shows 
the maximal obtained enhancement, the orange the mean enhancement of a generation respectively. 
The black dashed horizontal line marks the field enhancement of the reference dipolar antenna, 
b) shows the probability of each 21x21-matrix position for hosting a gold cube for the 10 best 
structures of the respective generation. The blue spot at the center marks the "gap", where the 
optimization took place and thus always was empty. 

features leads to a jump in the maximum fitness of a generation, when they are established 
the first time. The first such feature is a gap, which ensures the minimal possible distance 
of charge carriers to the optimization position. Then the areas next to the gap are getting 
more solid, to establish currents and supply the gap with a higher surface charge density. 
The first big jump in the fitness originates from the removal of the block directly above the 
center, opening the gap. 

The fitness also occasionally reduces which may occur, when a local maximum of the con- 
figuration space was found and among the offspring structures there is not yet an individual 
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that includes a new key feature. More and more individuals will then be generated by mu- 
tation, thus further lowering the fitness in the beginning but also generating new features, 
which will finally lead to the formation of new structural elements that will then generate 
better structures. 



QUALITY FACTOR 



The optimal antenna structure found by the evolutionary algorithm exhibits a compara- 
tively large quality factor Q ^ 20. Due to the eflFective wavelength scaling [2] the resonance 
length of plasmonic two- wire antennas strongly depends on the wire cross-section. A small 
cross section d leads to a reduced overall length for resonant antennas which in turn leads 
to a reduced radiation resistance. The reduced radiation resistance directly results in an 
increased quality factor of the antenna resonance. Table I demonstrates this behaviour for 
a series of antennas with the same resonance position, but varying cross section d. 




d / nm 


30 


25 


20 


15 


1 / nm 


90 


84 


74 


61 


X^/ nm 


725 


729 


729 


726 


AX 1 nm 


58 


46 


36 


29 


Q-factor 


12.5 


15.8 


20.3 


25.0 


F 

-^max 


-2200 


-3100 


-4100 


-4600 



TABLE L Q-factor, resonance width AA, and maximum near-field intensity enhancement £^max of 
optical antennas consisting of gold wires with length / separated by a 10 nm gap with constant 
resonance wavelength Ar but decreasing square cross section . 



Along with the higher Q-f actors, thinner antennas also show higher maximum near-field 
intensity enhancement. This is due to the better field confinement in the gap due to the 
smaller cross section and thus the enhanced lightning rod effect. 



RADIATION PATTERN 



To demonstrate the validity of the reciprocity theorem for the fittest antenna found by 
the evolutionary algorithm, we place a dipole emitter at the position of highest near-field 
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intensity enhancement and study its enhanced emission. The resulting far-field power fiow 
is depicted in Fig. 6. 



a) 




FIG. 6. Radiation patterns of reference and fittest matrix antenna, a) Definition of spherical 
coordinates 9 and Lp with respect to the antenna orientation. A point dipole emitter emitting at 
the resonance wavelength Ares — 650 nm is placed at the position of highest near-field intensity 
enhancement. Since the structures are symmetrical in z-direction, only the radiation power in 
the upper half space needs to be evaluated. The coloured markings (red, blue, green) denote 
particular constant-coordinate lines, b) shows geometry and normalized dipolar radiation pattern 
in cylindrical projection for the reference antenna, c) shows geometry and radiation pattern of the 
fittest matrix antenna normalized to the reference antenna. 
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The directivity of an antenna is defined as the maximum of the directive gain 

= (8) 

-'rad 

with Prad the integrated radiated power computed from the integral over the radiation power 
per unit sohd angle p{6^(p). The reference antenna shows a symmetrical dipolar radiation 
pattern with a directivity of Dref = 2.09 dBi which is only slightly smaller than Da/2 = 2.15 
dBi of an ideal thin wire half wave dipole. 

The radiation pattern of the best matrix antenna deviates only slightly from reference 
antenna, but exhibits a more pronounced radiation in the z-direction. Considering reci- 
procity, its enhanced directivity of Dea = 2.31 dBi reflects the optimization to irradiation 
by a gaussian beam. 
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